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We study the cos (fu azimuthal asymmetry in double polarized semi-inclusive pion production by 
considering dynamical twist-3 effects. In particular, we evaluate the role of the transverse momentum 
dependent distributions eL{x,k^) and gj^(x,k%) on the asymmetry. Using two different sets of 
spectator model results for these distributions, we predict the cosifh asymmetry of 7r~, and 7r° 
at the kinematic configuration available at CLAS, HERMES and COMPASS. Our estimates show 
that the asymmetries are positive for all the pions and could be accessed by CLAS and HERMES. We 
also find that gives the dominant contribution to the cos0/i asymmetry, while the contribution 
ol Cl is almost negligible. 

PACS numbers: 12.39.-x, 13.60.-r, 13.88.+e 


I. INTRODUCTION 

Nowadays it is very clear that a better understanding 
on the nucleon structure can be achieved if one goes be¬ 
yond the collinear picture to take into account the trans¬ 
verse degree of freedom of partons. Early investigation [ij 
on the unpolarized semi-inclusive deep inelastic scatter¬ 
ing (SIDIS) demonstrated that the intrinsic transverse 
motion of quarks can give rise to a cos (j)h asymmetric 
distribution of the final state hadron. This mechanism, 
usually called as the Cahn effect, provides a useful tool 
to probe the partonic transverse momentum which is still 
less understood so far, although the effect appears kine¬ 
matically at the subleading order of inverse hard scale. 
The same idea was then explored further by several ex¬ 
perimental and theoretical studies P-fl^ to give con¬ 
straints on the transverse structure of nucleon, i.e., the 
effect was applied [isini to extract the average values 
of the intrinsic transverse momenta of quarks inside the 
nucleon from SIDIS or Drell-Yan data [l^. 

Recently, in a new study H the Cahn effect was 
extended to the case of double longitudinally polarized 
SIDIS, prediciting a similar cos (j)h azmithual asymme¬ 
try that originates from the transverse momentum de¬ 
pendent (TMD) helicity distribution function giL^x, k^). 
Due to the positive value of and u quark dom¬ 

inance, the asymmetry based on the Cahn effect was 
found to be negative for both the charged and neutral 
pions in the case of proton target. 

In this work, we will study the cos (fh azimuthal asym¬ 
metry in double longitudinally polarized SIDIS in an al¬ 
ternative approach, that is, to employ dynamical twist-3 
effects. As shown in Ref. [l^, the polarized structure 


*Electronic address: zhunlu@seu.edu.cn 
^Electronic address: mabq@pku.edu.cn 


function that associated with the cos (f>h asymme¬ 

try can be expressed in terms of the twist-3 TMD distri¬ 
bution/fragmentation function combined with the twist- 
2 fragmentation/distribution function. Particularly, two 
twist-3 TMD distributions appear in the convolutions: 
the T-even distribution gj^{x,k^) and the T-odd distri¬ 
bution eij(x,k‘^). The former one can be decomposed 
into the following form via the equation of motion rela¬ 
tion [I^ 

xgL = xgL + aiL +j^Kl- (i) 

Taking the component giL from x gj; in the above equa¬ 
tion is equivalent to adopting the Cahn effect. In this 
study we will consider the effect of the entire twist-3 dis¬ 
tribution gj^. In addition, we also take into account the 
contribution of coupled with the Collins fragmenta¬ 
tion function [I3- We calculate these two TMD distri¬ 
butions of valence quarks inside the proton by employing 
the spectator diquark model and predict the correspond¬ 
ing cos (ph asymmetry for charged and neutral pions at 
the kinematics of JLab, HERMES and COMPASS. We 
note that sizable dynamical twist-3 effects may also ap¬ 
pear in other processes which involve different polariza¬ 
tions of the lepton beam and the nucleon target (l8l - [30l| . 


II. MODEL CALCULATIONS ON THE TWIST-3 
TMD DISTRIBUTIONS gp AND cl 

In this section, we briefly present our calculation on 
the distributions pp(x, k'^) and eL(x, using the spec¬ 
tator model [3ll - l^ We will consider the contributions 
from both the scalar diquark and the vector diquark. In 
the case of vector diquark components, we use two ap¬ 
proaches for comparison. The main differences between 
them are the form for the vector diquark propagator, as 
well as the flavor separation for u and d valence quarks. 



2 


The gauge-invariant quark-quark correlator for a lon¬ 
gitudinally polarized nucleon in SIDIS reads: 

Hx^kT) = I 

X /:€-[Or,|T]/:«"[oo-,r]V'*(C)|i"^L) • (2) 

Here the light-cone coordinate a = [a~,a'^,aT\ is em¬ 
ployed, £s are the gauge links ensuring the gauge invari¬ 
ance of the operator, and k and P are the momenta of the 
struck quark and the target nucleon, respectively. The 
distributions and thus can be obtained from the 
correlator m using the traces [H, [s^ 

-iTr[$7“75], (3) 

i Tr[$i75]. (4) 


Sl ^T^^P gj;{x,k^) = - 

„ M , 


Substituting © into o and © into ®, we obtain 
the following expressions for gj^ and Bl from the scalar 
diquark component: 


gi'^{x,kT) = - 


N^{1 - x)^ (1 - x)^M^ - M2 - 




{krp + 1/2)4 


(9) 


el(x,k‘^) = Cpa, 


Ns\l-xf (x + f)(L2-4) 

32^3 Ll{Ll + klf 


( 10 ) 


We note that the result in dH) has already been given in 

Ref. [m. 

The correlator contributed by the vector diquark can 
be obtained in the similar way which was applied to cal¬ 
culate and . Here we cast the expressions for 
the correlator from the vector diquark component at the 
lowest order: 


First we consider the contribution from the scalar di¬ 
quark. In the lowest order, the correlator may be cal¬ 
culated by suppressing the gauge link in the operator. 
After some algebra, we arrive at the expression for the 
correlator: 


^( 0 )/ , ^ + + M){lf: + m)\ 

* 327r3p+ (4+ £2)4 

(5) 

Here Ns is the normalization constant, and the notation 
£2 has the form 

£s = (1 — a:)A2 + xM"^ — x(l — x)M2, (6) 

with As being the cutoff parameter for the quark momen¬ 
tum and Ms the scalar diquark mass. 

The lowest order result in Eq. © can be used to cal¬ 
culate T-even distributions. However, it leads to a van¬ 
ishing result for T-odd distributions. In order to yield 
a nonzero contribution from the correlator, one has to 
consider the effect of the gauge links [H, Hj, Ii3| , or the 
rescattering between the struck quark and the specta¬ 
tor diquark. Here we expand the gange-links to the first 
nontrivial order, which corresponds to the one gluon ex¬ 
change approximation. At this order, the correlator has 
the form: 


$W(x,fcT) = igNl 


(I-x)2 




647r3(P+)2 (fc2 +L2)2 


d?qT [{li - i + + M){ + to)] 


(27r)2 


q2 [(kj, - + +2] 


(7) 


with = 0, and T^ being the vertex between the gluon 
and the scalar diquark: 


= ig{2P -2k + qY ■ (8) 


4>W( 


X, kr) = .. „ —a 


64^3p+ 


[{^ + + m)] 

{kr + 


( 11 ) 


and at the one-loop level: 

fer) = igNl J'l 


I287r3(P+)2 {kl + LlY 

X dcr/3{P - k + q) 

, [( ^ - ^ + "i)7'^75-?(Ar - If)Y{ ^ + ™)] 


Qt [(kr - QtY + Ll] 


( 12 ) 


respectively. In Eqs. (HU and (HU), we have used 
to denote the propagator of the vector diquark, which 
corresponds to the sum of its polarization vectors. Also, 
rM,a/3 (Jej^otes the vertex between the gluon and the vec¬ 
tor diquark 

= -ig[{2P -2k + qYg‘"^ - (P - fc + qYg^^^ 

- (P - fc) V“] ■ (13) 


In literature, different choices have been made for . 
As shown in Ref. [s^, different form of generally 
leads to different result of the correlator. In this work, 
we will consider two choices for for comparison. The 
first one has the form: 


d>^Yk) = -g^'" 


k^NL + k'^n^ 

k ■ n_ 


m2 

[fc ■ n_]^ 


(14) 


which is motivated by the light-cone formalism [s^ for 
the vector diquarks. Applying the propagator (fTH) . we 
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and we denote them as the Set I results of /^. 

The second form for the vector diquark propagator em¬ 
ployed in our calculation is 


d^''{k)= (17) 

which has been applied in Ref. . Similarly, using (EZD 
we obtain alternative expressions for and e\'. 


e^(a;, krp) 
g'^^{x, krp) 


N,^{l-x)^ix+§)iLl-k^p) 

Set II 32713 LULl + k^^)^ ’ 

(18) 


N^il - x)^ (1 - x)‘^M‘^ - m 2 - fe2 


Set II 


167r3 


{kp + 1/2)4 


(19) 


which we denote as Set II results. 

In order to obtain the flavor dependence of the TMD 
distributions, one should assign the relation between the 
quark flavors and the diquark types. In Ref. [s^l, a gen¬ 
eral relation is introduced: 

= Us., I = , (20) 

where a and a' represent the vector isoscalar diquark 
a{ud) and the vector isovector diquark a{uu), respec¬ 
tively, and Cs, Ca and Ca' are the parameters of the model. 
These parameters as well as the mass parameters (such 
as the diquark masses Mgjy, cut-off parameters Ag/y) 
are fitted from the ZEUS unpolarized parton distribu¬ 
tion functions 1411 an d GRSVOl polarized parton distri¬ 
bution functions |42j|. We combine Eqs. ®, (0, (HSl), 
and (HU) to obtain the Set I distributions and g^'^ 
with q = u and d. For the strong coupling appearing 
in the expressions for e^, we choose ~ 0.3. In the 
left panels of Figs. [T] and [U we plot the cc-dependence (at 
kr = 0.2 GeV) of the functions gp‘^{x, kp) and e|^(x, kp) 
timed with x for g = it and d quarks in Set I. We also 
plot the At- dependence (at x = 0.3) of the distributions 
in the right panels of Figs. [T]and[2J 

Different from Eq. dsni), another kind of flavor separa¬ 
tion has been employed previously [3l|, : 

/Is..,, = 5 /’+ 5 /“. /Is.. „ = /“'. ( 21 ) 

where the coefficients 3/2, 1/2 and 1 in front of are 
obtained from the SU(4) spin-flavor symmetry of the pro¬ 
ton wave function. In this model, the mass parameters 


FIG. 1: Model results for xgp{x, kp) in Set I. Left panel: the 
X dependence of xg^^ (solid line) and xgp‘^ (dashed line) at 
kp = 0.2 GeV; right panel: the kp dependence of xgp'^ (solid 
line) and xgp‘^ (dashed line) at x = 0.3. 
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FIG. 2: Similar to Fig. [2l but for the model results of xe/ 
(solid line) and xep (dashed line). 


for different types of vector diquarks are the same, and we 
apply the values for the parameters from Ref. |35l |. Us¬ 
ing the relation (ED, together with the expressions ([5]), 
(nni), (HU) and m , we obtain another set of TMD dis¬ 
tributions, which we denote as Set II distributions. The 
corresponding numerical results are plotted in Figs.[3]and 
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Comparing Fig. [T] with Fig. |3] and Fig. [2] with Fig. |3 
we can see that the sizes of the TMD distributions in 
Set I are different from those in Set II. In both sets, 
the signs of ez, and g^ turn to be negative in the spec¬ 
ified kinematics (x = 0.3 and kp = 0.2 GeV). Also, the 
size of the T-even distribution g^ is generally larger than 
that of the T-odd distribution cl- This is understand¬ 
able since T-odd distributions are yielded from the higher 
order expansion of gauge link. The distribution g^ has 


1.6 

. Set II 

1 

■ Set II 


/ \ k.p=0.2GeV 



0.8 


0 


0.0 


-1 

f' 


N. '' / 


x=0.3 

0.8 

'' - -X g vs. X 

-2 

■ / -X g / " vs. k.,. 

1.6 

- 

-3 

/ -X g vs. k.,- 


0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 


X MGeV) 

FIG. 3: Model results for xgp{x,kp) in Set II. 
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FIG. 4: Model results for xeL{x, k^) in Set II. 


also been calculated in the Bag model. We find that 
the tendency of the x dependence o f qj- in our calcu¬ 
lation agrees with the result in Ref. |43| . that is, there 
is a node in the intermediate x region. This behavior 
may be explained by the so called Lorentz invariance re¬ 
lation [l^, m, 115 . We also point out that our result 
for e\^{x, fcy) agrees with the time reversal constraint for 
distributions / (Pkx e1{x,k'^) = 0. 


III. PREDICTIONS ON THE cos 
ASYMMETRY FOR CHARGED AND NEUTRAL 
PIONS IN POLARIZED SIDIS 


In this section, we perform phenomenological analysis 
on the cos (f>h asymmetry for pions in SIDIS: 

r (£) + p-^{P) ^ + h{Ph) + X{Px ). (22) 

Here the arrow —>■ denotes the longitudinally polarization 
of the beam or proton target, (. and I' stand for the mo¬ 
menta of the incoming and outgoing leptons, and P and 
Ph denote the momenta of the target nucleon and the 
final-state hadron, respectively. The kinematics of SIDIS 
can be expressed by the following invariant variables 

P-q P-Ph 2Mx 

^ = 2p-q' y=T-v " = ^ = 

Q^ = -q\ s=(P + £)2, (23) 

with q = £ — £' the momentum of the virtual photon, and 
W the invariant mass of the hadronic final state. The 
reference frame adopted in this work is shown in Fig. [Sj 
in which the momentum of the virtual photon is along the 
z axis, and the longitudinal polarization of the target is 
along the opposite direction of 2 ; axis. Thus, we will not 
consider the contribution from the transverse component 
of the polarization, which involves the TMD distribution 
giT- In this frame, the transverse momentum of the final 
hadron with respect to the fragmenting quark is denoted 
by i-’r, and the azimuthal angle of the hadron around 
the virtual photon is defined as (ph- 

The differential cross section of SIDIS by scattering a 
longitudinally polarized lepton beam off a longitudinally 



FIG. 5: The kinematical configuration for longitudinally po¬ 
larized SIDIS process. The initial and scattered leptonic mo¬ 
menta define the lepton plane {x — z plane), while the momen¬ 
tum of the detected hadron together with the 2 : axis identifies 
the hadron production plane; the longitudinal spin of the nu¬ 
cleon is along the —2 axis. 


polarized target can be expressed as [l5l| 
do- ^ 27ra2 

dxdy dzdP!^dph xyQ'^ 2(1 — e)\ 2x) ^ 

-h S\\\e\/2l{l^^cosph H- I , 

(24) 


here Fuu and are the spin-averaged and spin- 

dependent structure functions, respectively, and the ra¬ 
tio of the longitudinal and transverse photon flux is de- 
fined by e = /4 ■ The ellipsis stands for the 

leading-twist double-spin asymmetry All which 

will not be analyzed in this work. 

By exploiting the notation 


C[wfD]=x'^^e^^ / d^kr / d^pTS^{zkT — Pt + Pt) 


X w{kT,PT)r{x,k'^)D‘^{z,p'^), 


(25) 


we express the structure functions Fjju and as [T^ 

Fuu=C[hD^], (26) 


f: 


iCOS (f>h 


LL 


2M 

~Q 


c 


-h 


Pt ■ h'l' 

M 


Pt ■ Pt 

zMh 


{xclH^) 


(27) 


Here, we introduce the unit vector Pt = and use Mh 
denote the mass of the final hadron. As we restrict our 
scope on the role of the twist-3 distributions in 
in Eq. (1271) we have suppressed the terms containing the 
twist-3 fragmentation functions E and . 

The asymmetry as a function of x can be cast 

into 

^ / dFj Jdy^dz Cf v^2g(l - e) Fl°l 

^ ’ J dP^ J dy J dz Cf Fuu 

( 28 ) 


























5 



#0.15 
8 =1 

0.10 

0.05 

0.00 

-0.05 

- 0.10 




t I I I I- . I . I. I-. 

0.2 0.4 0.5 0.6 0.3 0.6 0.9 

X z P,(GeV) 


15 

10 

05 

00 

■ 0 
- n 



srrr-. . 

— 

- 

05 


. • - Set 1 e^H\ 

• • • -Setlle^H; 



■ — -Set 1 

- Set [l 

10 


■ Set 1 total 

- Set II total 


- 1 . 1 

-.1.1. 

-.1.1.1 


0.2 0.4 

0.5 0.6 

0.3 0.6 0.9 


X 

Z 

P,(GeV) 


FIG. 6: Prediction on for (left panel), 7r“ (central panel), and tt® (right panel) vs x, z, and Pt in SIDIS at CLAS 

off a proton target. The thick and thin curves correspond to the asymmetries calculated from the TMD distributions in Set I 
and Set II, respectively. The dotted and dashed curves represent the asymmetries from the clH^ term and the gj^Di term, 
and the solid curves correspond to the total contribution. 


where the kinematical factor C/ is defined as 

In a similar way, we can define the z-dependent and the 
P^-dependent asymmetries. 

In order to estimate the numerical results of , we 

apply the model results of and e l obtained in the pre¬ 
vious section . As for the Collins function appearing 
in Eq. (Ell), we adopt the parametrization (the standard 
set) from Ref. [i^. To obtain the Collins function for 
neutral pion, we employ the following isospin relation: 

^ = 1 , ( 30 ) 

For the TMD fragmentation function Df which 

couples with the distribution gj^, we assume its pr de¬ 
pendence has a Gaussian form 

(31) 

t^\Pt) 

where D\{z) is the integrated fragmentation function 
D\{z), for which we adopt the leading order set of the 
DSS parametrization [d^. For the Gaussian width (p^), 
we choose its numerical value as 0.2 GeV^, following the 
fitted result in Ref. 0- Finally, in our calculation, 
we consider the kinematical constraints on the intrin¬ 
sic transverse momentum of the initial quarks given in 
Ref. [13. 

We apply the following kinematics to estimate the 
asymmetry at CLAS: 

Ee = 5.5 GeV, 0.1 < cc < 0.6, 0.4 < z < 0.7, 

> 1 GeV^ Pt > 0.05 GeV, > 4GeV^. 

In the left, central, and right panels of Fig.[6l we plot the 
asymmetry for 7r+, 7r“ and 7r° as functions of ai, z, and 
Pt- The thick and thin curves correspond to the asym¬ 
metries that calculated from the TMD distributions of 


Set I and Set II, respectively. The dotted curves show 
the asymmetries contributed by e^, the dashed curves 
show those contributed by gj;, while the solid curves de¬ 
note the total contribution. We find that the asymme¬ 
tries calculated from both sets of TMD distributions are 
positive for all three pions. Also, It is clear that the asym¬ 
metries contributed by the T-even distribution gj^ domi¬ 
nate, and those contributed by are almost negligible. 
The asymmetries calculated from Set I TMD distribu¬ 
tions are about 5 to 10 percent in magnitude, and are 
several times larger than those from Set II distributions. 

Using the two sets of TMD distributions, we also pre¬ 
dict the asymmetry at HERMES with a 27.6 GeV 

positron beam off a proton target (Tsj l: 

0.023 <x< 0.4, 0.1 < y < 0.85, 0.2 < z < 0.7 

> lOGeV^ > lGeV^ 4 GeV < Ph < 13.8 GeV, 

and at COMPASS with a 160 GeV muon beam scattered 
off a deuteron target pol| - 

0.004 < a; < 0.7, y > 0.1, 0.2 < z < 0.9, 
a;F > 0, > 1 GeV^ O.lGeV < Pt < 1 GeV, 

5 GeV < W < 18 GeV. 

In Figs. [7] and [H we show our prediction on the cos (ph 
asymmetry for charged and neutral pions at HERMES 
and COMPASS, respectively. We find that the asymme¬ 
try at HERMES is smaller than that at CLAS. Again, 
the T-even distribution gj^ gives the dominant contribu¬ 
tion. The asymmetries for all pions from the two sets 
of TMD distributions at COMPASS are consistent with 
zero (less than 0.5%). 

IV. CONCLUSION 

In this work, we investigated the cos (ph azimuthal 
asymmetry in the double longitudinally polarized SIDIS. 
Particularly, we focused on the role of the genuine twist- 
3 TMD distributions and ignored the contribution from 
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FIG. 7: 
target. 


Similar to Fig. |6] but for of tt'*', tt , and 7r° as functions of x, z and Pt in SIDIS at HERMES off a proton 



FIG. 8: The asymmetry for tt'*', tt , and tt® as functions of x, z and Pt in SIDIS at COMPASS off a deuteron target. 


the twist-3 fragmentation functions. To give a quantita¬ 
tive estimate on the cos (ph asymmetry for different pions, 
we calculated e l and of the valence quarks within the 
framework of spectator model. We considered two differ¬ 
ent forms for the propagator of the vector diquark as well 
as different choices on the flavor separation to obtain two 
sets of TMD distributions, which were used to predict the 
asymmetry for 7r+, tt~ and at the kinematics 

of CLAS, HERMES and COMPASS. We found that the 
asymmetries from both sets of TMD distributions are 
sizable at CLAS, while at HERMES only the asymmetry 
from the distribution of Set I is measurable. Therefore, it 
would be feasible to access the cos ph asymmetry at least 
at CLAS. We also found that the predicted asymmetries 
are positive for all the pions. This is different from the 
estimate [ij based on the Cahn effect, which predicts 
negative asymmetries for pions, although the magnitude 
of the asymmetry in our calculation is consistent with 


the result in Ref. M- A positive cos cph asymmetry thus 
can be viewed as a clear signal of dynamical twist-3 effect 
that is different from the Cahn effect. Furthermore, our 
study shows that the main contribution to the asymme¬ 
try is from the T-even distribution gj^, while the contri¬ 
bution from the T-odd distribution cl almost vanishes. 
Future experimental data on the azimuthal asymmetry 
in the double longitudinally polarized reaction will clar¬ 
ify the role of the twist-3 TMD distributions. 
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